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Abstract 

A  multiple-slip  dislocation-density-based  formulation  and  computational 
schemes  that  are  coupled  to  grain-boundary  (GB)  interfacial  schemes  and  an 
internal  porosity  formulation  are  used  to  analyse  the  behaviour  and  interaction 
of  different  arrangements  and  geometries  of  explicit  pairs  of  voids  in  a 
polycrystalline  fee  aggregate.  The  GB  regions  are  treated  as  regions  with 
properties  and  topologies  that  are  distinct  from  that  of  the  grain  bulk.  The  GB 
kinematic  scheme  accounts  for  dislocation  density  interactions  with  GBs,  such  as 
dislocation  density  impedance,  blockage  and  GB  absorption.  These  evolving 
interfacial  conditions  are  monitored  throughout  the  deformation  history.  The 
analysis  indicated  that  void-to-void  interactions  result  in  dislocation  density 
evolution  and  saturation  and  porosity  localization  that  are  intricately  related  to 
both  dislocation  density  pile-ups  and  blockages  at  GB  interfaces,  and  GB 
absorption  within  different  GB  regions. 


§1.  Introduction 

GB  structure  and  orientation  are  crucial  factors  in  the  characterization  of 
ductile  and  brittle  failure  modes  in  polycrystalline  materials.  These  GBs,  generally 
consisting  of  transition  regions  of  highly  distorted  crystalline  structure  between  mis- 
oriented  grains,  with  material  properties  that  are  generally  distinct  from  those  of 
bulk  grains,  can  accelerate,  retard  or  inhibit  the  initiation  and  the  evolution  of  failure 
modes.  As  noted  by  Schmitz  et  al.  (1989)  and  Davies  and  Randle  (2001),  the  early 
stages  of  material  failure  in  crystalline  materials  are  not  characterized  by  crack 
nucleation  in  individual  grains,  but  the  critical  first  step  is  the  propagation  of  the 
crack  across  or  along  a  GB  plane.  The  dislocation  characteristics,  the  interface  and 
grain  properties,  the  topological  effects  of  internal  microstructure,  and  the  cohesive 
strength  of  evolving  internal  GB  properties  are  intricately  related  to  failure  initiation 
and  evolution. 

The  physical  scale  on  which  GB  effects  are  investigated  has  a  direct  bearing 
on  understanding  and  accurately  predicting  the  material  mechanisms  that  trigger 
failure.  On  a  macroscopic  spatial  scale,  it  is  generally  assumed  in  most  analyses  of 
polycrystalline  deformations  that  the  grains  are  homogeneous,  and  that  the  GBs 
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are  only  locations  where  stress  equilibrium  and  strain  boundary  conditions  have  to  be 
satisfied.  The  presence  of  GBs,  on  this  scale,  results  in  the  redistribution  of  stresses 
due  to  compatibility  requirements  for  elastic  and  inelastic  deformations.  However, 
since  crystalline  grains  are  generally  anisotropic,  the  deformations  due  to  external 
stresses  are  not  generally  compatible,  and  therefore  additional  stresses  must 
be  induced  to  match  the  shape  of  a  strained  grain  with  neighbouring  grains. 
Furthermore,  when  the  grains  are  plastically  deformed  by  inelastic  slip  on  specific 
crystallographic  slip  systems,  the  deformations  for  dilferent  orientations  are  not 
compatible  and  additional  stresses  may  also  arise.  Plastic  slip  is  also  generally  inho¬ 
mogeneous,  since  it  is  spatially  localized;  this  can  further  redistribute  the  stresses 
at  the  GB  (for  example  Kroner  (1986)).  On  a  microscopic  scale,  lattice  dislocations 
interact  with  the  GB,  which  acts  as  an  obstacle  to  their  motion.  The  strength  of 
these  obstacles  depends  on  the  interaction  of  the  GBs  with  the  dislocations.  The 
dislocations  impinging  on  the  GB  can  increase  the  magnitude  of  the  stress  field  in 
the  adjoining  grains  and  may  promote  nucleation  of  lattice  dislocations  in  the 
grain.  When  dislocations  penetrate  from  one  grain  to  the  next,  residual  dislocations 
may  be  left  behind  in  the  boundary.  Therefore,  the  strength  of  the  obstacle  may 
depend  on  the  Burgers  vector  of  residual  boundary  dislocations  that  are  created  by 
penetrating  lattice  dislocations  (Gleiter  1982). 

Furthermore,  experimental  results  (for  example  Lee  et  al.  (1990)  and  Marguiles 
et  al.  (2001))  have  clearly  indicated  how  slip  transmits  through  GBs  and  how,  if  the 
slip  is  impeded,  cracks  can  nucleate  along  the  GB  plane,  resulting  in  intergranular 
fracture.  The  significance  of  these  results  is  that,  owing  to  plastic  deformation,  dislo¬ 
cations  pile  up  against  an  interface,  and  then  cracks  nucleate  and  grow  along 
the  GB.  Therefore,  it  is  essential  not  only  to  deal  with  a  density  of  dislocations 
for  large  inelastic  deformations  in  aggregates  with  different  GB  interfaces  but  also 
to  treat  GBs  as  interfacial  planes  with  properties  distinct  from  those  of  the  bulk, 
where  a  population  of  dislocations  from  different  adjacent  grains  are  absorbed  or 
nucleated. 

One  of  the  current  challenges  is  how  to  account  for  different  dislocation 
density  interactions  with  GB  interfaces  with  distinct  topologies  and  properties, 
and  how  these  interactions  affect  ductile  failure  nucleation  and  growth  on  different 
physical  scales.  Interfacial  misorientations  between  grains  due  to  the  translation 
and  rotation  of  one  crystalline  lattice  with  respect  to  adjacent  grains  in  combination 
with  the  other  defects  such  as  cracks  and  voids  can  result  in  the  formation  of  dislo¬ 
cation  density  pile-ups  at  GB  interfaces,  or  the  absorption  of  dislocation  densities 
within  the  GB  region,  or  the  partial  transmission  of  dislocation  densities  across 
GB  interfaces  to  neighbouring  grains  (Watanabe  1989,  Baker  and  Liu  1994). 
Numerous  investigations  have  clearly  indicated  that  intrinsic  thermomechanical 
properties,  operative  crystallographic  systems  and  dislocation  mechanisms  asso¬ 
ciated  with  GB  regions  all  may  be  substantially  different  from  those  associated 
with  the  aggregate  crystalline  bulk  response  (Dingley  and  Pond  1979,  Baker  and 
Liu  1994,  Randle  1997,  Davies  and  Randle  2001).  However,  the  effects  of  how 
these  GB  interactions  and  mechanisms,  such  as  slip  transmission,  impedance  and 
blockage,  affect  void  nucleation  and  growth  on  different  physical  scales  are  not 
well  understood  and  have  not  been  adequately  quantified  (Watanabe  and 
Tsurekawa  1999). 

The  major  objective  of  this  study  is  to  understand  how  dislocation  density 
interactions  with  GB  interfaces  affect  microvoid  nucleation  and  growth,  and  how 
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this  is  related  to  material  behaviour  and  interactions  with  different  explicit  pairs  of 
void  distributions  and  arrangements  at  the  local  and  global  levels.  A  kinematically 
based  scheme  that  is  coupled  to  a  dislocation  density  multiple-slip  formulation 
(Ashmawi  and  Zikry  2002)  is  used  to  account  for  the  interactions  of  dislocation 
densities  with  GB  interfaces  in  polycrystalline  aggregates.  This  methodology 
is  used  with  an  internal  porosity  formulation  that  accounts  for  microvoid 
growth  and  nucleation.  In  this  study,  the  microvoid  nucleation  internal  variable 
is  a  function  of  updated  values  of  the  total  dislocation  density,  stress  triaxiality, 
accumulated  plastic  strains  and  temperature.  GB  interfaces  are  treated  as  regions 
with  distinct  width  properties  and  topologies.  Schemes  are  then  introduced  to 
monitor  how  dislocation  density  evolutions  on  different  operative  slip  systems 
interact  with  GB  interfaces,  and  how  dislocation  density  evolution  can  result  in 
slip  transmission,  impedance  and  blockage.  Furthermore,  it  is  shown  how  this 
can  result  in  dislocation  density  pile-ups  at  GB  interfaces,  or  dislocation  density 
absorption  within  the  GB  region,  or  dislocation  density  transmission  to  neighbour¬ 
ing  grains. 

This  paper  is  organized  as  follows:  the  dislocation-density-based  multiple-slip 
crystalline  formulation  is  introduced  in  §2;  the  dislocation  density  evolution  is 
presented  in  §3;  the  kinematic  scheme  for  the  interaction  of  dislocation  densities 
and  slip  with  GB  interfaces  is  outlined  in  §4;  the  internal  porosity  formulation  is 
introduced  in  §5;  the  computational  method  is  given  in  §6;  the  results  are  presented 
and  discussed  in  §7;  the  summary  and  conclusions  are  given  in  §8. 


§2.  Dislocation  density-based  multiple-slip  constitutive  formulation 
The  formulation  for  the  multiple-slip  crystal  plasticity  rate-dependent  constitu¬ 
tive  relations,  and  the  derivation  of  the  evolutionary  equations  for  the  mobile  and 
immobile  dislocation  densities,  which  are  coupled  to  the  multiple-slip  crystalline 
formulation,  are  presented.  The  formulation  is  based  on  that  developed  by  Zikry 
and  Kao  (1996)  and  Ashmawi  and  Zikry  (2000,  2003). 

It  is  assumed  that  the  velocity  gradient  can  be  decomposed  into  a  symmetric  part, 
the  deformation  rate  tensor  Dy  and  an  antisymmetric  part,  the  spin  tensor  Wy.  It  is 
further  assumed  that  the  total  deformation  rate  tensor  Dy  and  the  total  spin  tensor 
Wy  can  be  then  additively  decomposed  into  elastic  and  plastic  components: 

Dy  -  Dy  +  Z)P,  (la) 

Wy  —  Wy  +  Wy,  (lb) 

in  which  Wy  includes  the  rigid-body  spin.  The  inelastic  parts  are  defined  in  terms  of 
the  crystallographic  slip  rates  as 


dp  _  pf-y-),  (2fl) 

w*  -  cofY(a\  (2b) 

where  a  is  summed  over  all  slip  systems,  and  P{y]  and  off  are  second-order  tensors, 
defined  in  terms  of  the  unit  normals  to  the  slip  planes  and  the  unit  slip  vectors  to  the 
slip  directions. 
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For  rate-dependent  inelastic  materials,  the  constitutive  description  on  each  slip 
system  can  be  characterized  by  a  power-law  relation 


r(«) 

Aa)  _  Aft)  x_ 
y  ^ref  (a) 

Tref 


no  sum  on  a, 


(3) 


where  y(rff  is  the  reference  shear  strain  rate  which  corresponds  to  a  reference 
shear  stress  and  m  is  the  rate  sensitivity  parameter.  The  reference  stress  that  is 
used  here  is  a  modification  of  widely  used  classical  forms  (Mughrabi  1987)  that 
relate  the  reference  stress  to  a  square  root  dependence  on  the  dislocation  density  as 


^  =  ^  +  ab±al(^)'12,  ,4, 

where  G  is  the  shear  modulus,  b  is  the  magnitude  of  the  Burgers  vector,  tyff)  is  the 
static  yield  stress  and  the  coefficients  a j  are  interaction  coefficients  that  generally 
have  a  magnitude  of  unity. 


§  3.  The  evolutions  of  mobile  and  immobile  dislocation  densities 
At  a  given  state  for  a  deformed  material,  the  dislocation  structure  of  total 
dislocation  density  p<a)  can  be  assumed  to  be  additively  decomposed  into  a  mobile 
dislocation  density  p^*  and  an  immobile  dislocation  density 

p^pW  +  pM  (5) 


It  is  assumed  that,  during  an  increment  of  strain,  an  immobile  dislocation  density 
rate  is  generated  and  an  immobile  dislocation  density  rate  is  annihilated  for  statis¬ 
tically  stored  dislocation  densities  (for  example  Liu  et  al.  (1998)).  The  balance 
between  dislocation  generation  and  annihilation  equations  is  the  basis  for  the 
evolution  of  mobile  and  immobile  dislocation  densities  as  a  function  of  strain. 
Based  on  these  arguments,  it  can  be  shown  (see  the  paper  by  Kameda  and  Zikry 
(1996)  for  a  detailed  presentation)  that  the  coupled  set  of  nonlinear  evolutionary 
equations  of  mobile  and  immobile  dislocation  densities  can  then  be  given  by 


_  Aa) 

dt  y 


ksour 

IF 


Pim  km  in  ter  (  7/  \  kimmob  /  (a)\  *7“ 
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dt  y 
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^recov 


(6) 

(7) 


where  gsour  is  a  coefficient  pertaining  to  an  increase  in  the  mobile  dislocation  density 
due  to  dislocation  sources,  kminter  is  a  coefficient  related  to  the  trapping  of  mobile 
dislocations  due  to  forest  intersections,  cross-slip  around  obstacles  or  dislocation 
interactions,  igrccHV  is  a  coefficient  related  to  the  rearrangement  and  annihilation 
of  immobile  dislocations,  kimmob  is  a  coefficient  related  to  the  immobilization  of 
mobile  dislocations,  H  is  the  activation  enthalpy  and  k  is  Boltzmann’s  constant. 
As  these  evolutionary  equations  indicate,  the  dislocation  activities  related  to  recov¬ 
ery  and  trapping  are  coupled  to  thermal  activation. 

To  couple  the  evolutionary  equations  to  the  multiple-slip  crystal  plasticity 
formulation,  the  four  g  coefficients  in  equations  (6)  and  (7)  and  the  enthalpy  H 
must  be  determined  as  functions  of  the  deformation  mode.  The  enthalpy  H  is 
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determined  by  defining  an  exponential  ratio  of  the  current  temperature  to  the  refer¬ 
ence  temperature.  The  four  g  coefficients  are  determined  by  using  the  following 
two  general  conditions  pertinent  to  the  evolution  of  dislocation  densities  in  crystal¬ 
line  materials:  that  the  mobile  and  immobile  densities  saturate  at  large  strains, 
and  that  the  relaxation  of  the  mobile  dislocation  density  to  a  quasisteady-state 
value  occurs  much  more  rapidly  than  the  variation  in  the  immobile  density 
(Mughrabi  1987,  Hansen  1990,  Bay  et  al.  1992). 


§4.  Kinematics  of  dislocation  density  interactions  with 

GRAIN  BOUNDARIES 

As  stated  earlier,  GBs  can  act  as  effective  barriers  to  motion  of  lattice  disloca¬ 
tions.  Dislocation  densities  can  accumulate  within  GB  regions,  and  this  accumula¬ 
tion  within  GB  interfacial  regions  can  lead  to  substantial  increases  in  GB  internal 
stresses  and  GB  misorientations.  As  GB  misorientations  approach  large  angle 
values,  GBs  act  as  barriers  to  crystalline  slip  in  the  lattice,  which  can  result  in 
the  formation  of  pile-ups  that  subsequently  has  a  direct  consequence  on  failure 
initiation  and  evolution. 

There  is  a  myriad  of  dislocation  density  interactions  with  GB  interfaces,  such  as 
GB  absorption  of  lattice  dislocations  without  dissociation  into  grain-boundary 
dislocations  (GBDs)  (Baker  et  al.  1987,  Baker  and  Liu  1994),  partial  dislocation 
transmission  from  one  grain  to  the  adjacent  grain  with  a  residual  GBD  left  in  the 
GB  region  (Lee  et  al.  1992,  Baker  and  Liu  1994),  full  dislocation  transmission 
from  one  grain  to  the  adjacent  grain  with  no  residual  GBD  left  in  the  GB  region 
(Shen  et  al.  1988,  Lee  et  al.  1992,  Baker  and  Liu  1994),  and  dislocation  absorption 
and  subsequent  re-emission  from  the  GB  (Shen  et  al.  1988). 

In  this  study,  we  have  focused  on  investigating  the  following  representative 
interactions: 

(i)  full  and  partial  dislocation  density  transmission  from  one  grain  to  neigh¬ 
bouring  grains; 

(ii)  full  and  partial  dislocation  density  transmission  into  a  GB  and  blockage  at 
neighbouring  grains; 

(iii)  dislocation  density  impedance  and  potential  pile-ups. 

These  kinematic  scenarios  provide  a  general  methodology  that  can  be  used  as  a 
framework  for  physically  representative  GB  interfacial  mechanisms  and  interactions. 

A  detailed  presentation  of  the  kinematic  scheme  for  dislocation  density  trans¬ 
mission  and  impedance  from  one  grain  to  neighbouring  grains  across  GB  regions 
has  been  given  by  Ashmawi  and  Zikry  (2002).  Only  a  brief  outline  will  be  given  in 
this  paper.  The  following  steps  are  applied  to  determine  the  transmission,  impe¬ 
dance,  and  absorption  of  dislocation  densities  at  each  GB  interface  in  the  aggregate. 

(i)  A  potential  dislocation  density  envelope  is  used  to  identify  regions  of  high 
dislocation  density  activity  (figure  1)  on  each  corresponding  slip  system  in 
regions  adjacent  to  all  GB  interfaces.  This  dislocation  density  envelope 
encompasses  all  the  finite  elements  immediately  adjacent  to  the  GB  and 
that  extend  to  the  centre  of  each  grain. 

(ii)  Coordinate  frame  transformations  are  performed  with  respect  to  each 
adjacent  GB  interface  to  identify  the  direction  of  the  slip  system  orientation. 
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|  Leading  Dislocation  Density  Element 
J  Potential  Dislocation  Density  Pileup  Envelope 
^  GB  Potential  Transmission  Zone 

|  Neighbouring  Grain  Potential  Transmission  Zone 

Figure  1 .  Dislocation  density  envelopes;  lp  is  the  radial  pile-up  length. 


(iii)  Geometric  compatibility  conditions  are  then  obtained  to  determine  whether 
neighbouring  slip  systems  are  oriented  for  transmission  or  impedance 
through  the  GB. 

(iv)  If  transmission  is  then  possible  through  the  GB,  a  transformation  factor  is 
calculated  to  provide  a  measure  of  how  much  of  the  dislocation  density 
penetrates  through  the  GB  and  into  neighbouring  grains. 


4.1.  Potential  dislocation  density  envelope 
It  is  essential  to  determine  how  dislocation  densities  evolve  and  accumulate 
near  GB  regions  and  what  the  maximum  values  and  corresponding  active  slip 
system  directions  and  orientations  are,  such  that  slip  system  impedance  and  trans¬ 
mission  can  be  controlled  and  monitored  in  polycrystalline  aggregates.  Slip  impe¬ 
dance  at  a  GB  can  lead  to  pile-ups  and  stress  accumulation,  if  neighbouring 
slip  systems  within  the  GB  and  adjacent  grains  are  not  activated  and  oriented  for 
slip  transmission.  Therefore,  it  is  crucial  to  identify  the  critical  regions  of  dislocation 
density  accumulation.  A  potential  dislocation  density  envelope  is  used  to  identify 
these  critical  regions,  and  it  is  determined  as  follows. 

(i)  The  maximum  dislocation  density  on  the  corresponding  slip  system  is  iden¬ 
tified  in  the  region  adjacent  to  the  GB  by  monitoring  the  slip  rate  and  dislo¬ 
cation  density  activity  on  each  slip  system. 

(ii)  This  process  is  repeated  until  half  of  each  grain  diameter  is  transversed. 
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From  this  information,  the  envelope  is  constructed  as  shown  in  figure  1.  This 
radial  distribution  is  not  constrained  to  be  of  any  predetermined  mathematical  form. 
What  this  information  provides  is  a  distribution  of  the  total  dislocation  density  as  a 
function  of  length  orientation.  This  length  /p  is  a  radial  line  that  emanates  from  the 
GB  to  the  middle  of  the  grain.  This  radial  line  encompasses  the  dislocation  density 
envelope.  This  potential  dislocation  density  envelope  can  now  be  used  to  identify 
regions,  in  all  eight  adjacent  grains,  of  potential  transmission,  or  impedance,  or 
blockage. 


4.2.  Dislocation  density  orientation,  grain-boundary  coordinate  frame 
transformation  and  slip  plane  geometric  compatibility 
To  ascertain  how  dislocation  densities  evolve  near  GB  interfaces,  it  is  essential  to 
determine  whether  dislocation  densities  are  moving  away  or  towards  the  GBs. 
Flence,  the  orientation  and  direction  of  slip  with  respect  to  the  GB  will  have  to 
be  determined.  The  GB  normal  vector  iigb  and  the  in-plane  vector  mgb  are  used  to 
determine  dislocation  density  directions  and  orientations.  These  coordinate  transfor¬ 
mations  have  to  be  performed  for  all  slip  systems  on  all  eight  GB  interfaces  for  each 
grain. 

Dislocation  density  transmission  or  impedance  across  GB  interfacial  regions 
are  dependent  not  only  on  whether  there  are  adjacent  active  slip  systems  but 
also  on  whether  these  slip  systems  are  geometrically  compatible,  such  that  disloca¬ 
tion  density  transmission  is  possible  across  different  slip  systems.  In  this  analysis, 
geometric  compatibility  will  be  defined  by  two  angles  6  and  ft  (figure  2).  If  these 
values  are  exceeded,  then  transmission  through  the  GB  is  not  possible.  The  angle 


Figure  2.  General  slip  plane  geometric  compatibility. 
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between  the  two  intersection  lines  (vectors,  1]  and  12)  that  represent  the  lines 
formed  by  the  intersection  of  slip  planes  and  the  GB  interfacial  layer  is  given  by 


cos  0'  — 


li  12 

|li  I  II2I 


(8) 


and  0  is  determined  on  the  basis  of  the  following  two  conditions: 


0  = 


W\ 

if 


(9) 


Geometrically,  it  is  clear  that  slip  transmission  is  a  maximum  when  0  =  0°.  In 
this  case,  full  transmission  would  occur.  However,  as  noted  by  Werner  and  Prantl 
(1990)  and  Davies  and  Randle  (2001),  as  0  increases,  dislocation  density  transmission 
will  be  impeded.  Also,  slip  planes  intersecting  the  boundary  generally  maintain  an 
angle  of  15°  or  less  (Werner  and  Prantl  1990).  The  second  condition  for  ensuring 
geometric  compatibility  is  therefore  to  monitor  both  slip  plane  orientations  along  the 
interfacial  region  with  respect  to  an  axis  normal  to  a  plane  that  includes  GB  normals 
and  in-plane  vectors.  An  angle  ft  for  this  geometric  compatibility  can  be  determined 
by  defining  the  dot  product  of  two  normals  fi]  and  n2  to  the  slip  planes  on  either  side 
of  the  GB  plane  as 

cos  ft  —  fij  •  n2.  (10) 


Full  transmission  (with  respect  to  ft  compatibility)  of  slip  takes  place  when  ft  =  0. 
As  a  kink  starts  to  form,  for  increasing  value  of  0  =  0°,  dislocation  densities  will 
be  impeded.  We  shall  assume  that  ft  compatibility  is  bounded  by  a  value  of  35°. 
As  noted  by  Werner  and  Prantl  (1990),  these  critical  values  of  slip  orientation 
are  based  on  experimentally  observed  range  of  values  for  slip  compatibility.  In 
summary,  dislocation  densities  can  transmit  when  conditions  of  both  0  compatibility 
and  ft  compatibility  are  satisfied  in  conjunction  with  the  other  kinematic  conditions 
outlined  in  this  section. 


4.3.  Transmission  distribution  factor 

Geometric  slip  system  compatibility  affects  to  a  larger  degree  dislocation  density 
transmission.  The  most  favourable  configuration  occurs  when  0  =  0°  and  p  —  0°, 
which  indicates  complete  transmission.  Deviation  from  this  configuration  would 
represent  a  partial  transmission  through  the  interface.  A  transmission  factor  £, 
based  on  0  and  p,  will  be  defined  as 

ft  —  cos  0  cos  p.  (11) 

This  will  provide  a  measure  of  how  much  of  the  dislocation  density  penetrates 
through  the  GB  and  into  neighbouring  grains.  All  activities  on  all  slip  systems  have 
to  be  continually  monitored  on  all  eight  GB  regions  for  each  grain  in  the  aggregate. 
In  some  cases,  pile-ups  can  be  relieved  by  slip  system  rotations  in  neighbouring 
grains,  or  pile-ups  may  intensify  as  the  deformation  and  failure  modes  evolve.  In 
other  cases,  dislocation  density  transmission  through  the  GB  may  occur  by  partial 
transmission  on  different  slip  systems  into  neighbouring  grains.  All  these  situations 
need  to  be  delineated  for  accurate  predictions  of  the  effects  of  GB  interfacial  regions 
on  overall  behaviour. 
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§  5.  Internal  void  porosity  evolution 

As  noted  earlier,  interrelated  physical  mechanisms  can  result  in  void  nucleation, 
growth  and  coalescence.  In  this  section,  an  internal  porosity  evolution  is  introduced. 
This  evolution  relation  will  be  based  on  an  internal  variable  /,  which  will  account  for 
microvoid  nucleation  and  growth.  This  internal  porosity  relation  will  be  coupled 
with  the  proposed  dislocation-density-based  crystalline  constitutive  formulation 
and  the  interfacial  GB  dislocation  density  interaction  kinematic  models  to  obtain 
predictions  of  void  nucleation,  growth  and  coalescence  in  polycrystalline  aggregates. 
This  will  provide  an  understanding  of  how  dislocation  density  evolution  and  GB 
interfaces  affect  microvoid  nucleation  and  growth  and  how  this  is  related  to  interac¬ 
tions  with  discrete  void  distributions.  This  will  provide  failure  scenarios  that  can 
be  used  to  determine  how  explicit  voids  interact  with  GB  and  dislocation  density 
activities. 

Void  failure  modelling  will  encompass  two  concurrent  approaches.  One 
approach  will  be  based  on  the  internal  microvoid  variable  /,  which  will  be  coupled 
to  the  proposed  methodologies,  and  the  second  approach  will  be  based  on  modelling 
distributions  of  explicit  voids  with  different  geometries  and  morphologies  within  the 
aggregate. 

Interrelated  physical  mechanisms  can  result  in  void  nucleation,  growth  and 
coalescence.  The  internal  porosity  evolution  relation  is  based  on  an  internal  variable, 
which  accounts  for  microvoid  nucleation  and  growth.  It  is  assumed  that  the  porosity 
evolution  is  due  to  microvoid  growth  and  nucleation: 

f  ./growth  +  /nucleation  •  (12) 

The  rate  /growth,  of  the  growth  of  existing  microvoid  volume  fraction  is  assumed  to  be 
due  to  matrix  incompressibility: 

/growth  =  (1  ~f)Dlk,  (13) 

where  Dvkk  is  the  trace  of  the  inelastic  part  of  the  deformation  rate  tensor. 

We  assume  that  microvoid  nucleation  is  a  function  of  stress  triaxiality,  plastic 
slip  accumulation,  temperature  (Garrison  and  Moody  1987,  Magnusen  et  al.  1988, 
Needleman  1989,  Becker  and  Smelser  1994)  and  dislocation  density  activities 
(Garrison  and  Moody  1987).  Based  on  the  discussion  given  in  §4  and  detailed  in 
the  paper  by  Ashmawi  and  Zikry  (2002),  these  dislocation  density  activities  can 
encompass  mechanisms  such  as  slip  incompatibilities  and  intersections  that  can 
lead  to  stress  accumulation  and  pile-ups.  It  will  be  assumed  that  these  variables 
have  bounded  nucleation  rates  as  follows: 

* - >  h,  (14 a) 

Ptsat 


m  / 
-  — ►  '2> 


(14  b) 


Y  -*  h,  (14  c) 

^  ->  U,  (14 d) 
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where  p,  is  the  total  dislocation  density  for  the  most  active  slip  system  /l,  ptsat  is  the 
total  saturated  dislocation  density,  am  is  the  mean  normal  stress  (au/ 3),  ay  is  the 
static  yield  stress,  y  is  the  accumulated  plastic  shear  strain,  T  is  the  current  tempera¬ 
ture  and  Tre f  is  the  reference  temperature.  These  variables  will  be  combined  as 


A  = 


1 

/|  hhU 


(15) 


Since  nucleation  is  a  continuous  random  process,  the  evolution  formulation 
is  assumed  to  be  a  statistical  function  of  these  pertinent  parameters.  Therefore, 
the  microvoid  nucleation  rate  can  be  postulated  to  evolve  as 

/nucleation  =  *( b)v(am,  ,  T.  K)/nl,  (16) 

where  fni  is  a  nucleation  rate  limit,  and  the  random  variable  ¥  will  be  given  by 

08)  rr 

¥  =  A  y,  (17) 

Ptsat  ref 


such  that  0  ^  ¥  ^  1,  where  ¥  is  a  field  variable  that  represents  a  mapping  of  the 
interrelated  local  microstructural  quantities  that  affect  microvoid  nucleation  and 
growth.  Since  the  internal  porosity  is  generally  random,  it  will  be  assumed  that  ¥ 
is  a  continuous  random  variable  with  gamma  distribution,  such  that  the  distribution 
of  ¥  can  be  represented  with  a  probability  distribution  function  of 


X(¥;  a,  b )  = 


b‘T(a) 


i pa~ 


(18) 


where  T(a)  is  the  gamma  function  of  a.  The  parameters  a  and  b  have  positive  values. 

The  nucleation  rate  limit  /nl  is  defined  to  be  an  upper  material  limit  to  microvoid 
nucleation 


fn\  =  (fcr-f)ld[-  (19) 

where  /cr  is  the  maximum  porosity  for  a  material,  f  is  the  updated  porosity  and  d t 
is  the  incremental  time  step. 

The  updated  porosity  is  coupled  with  the  flow  stress  associated  with  each  slip 
system  a  as 

tJ  =  ^fd  “A  (20) 

where  f  ^  is  the  effective  reference  shear  stress.  Parametric  studies  were  performed 
to  assess  the  effects  of  varying  these  parameters  and  the  initial  porosity  together  with 
the  variables  a  and  b.  Limiting  values  for  each  of  the  variables  are  given  in  table  1. 


§  6.  Computational  techniques 

The  total  deformation  rate  tensor  Dtj  and  the  plastic  deformation  rate  tensor  D? 
are  needed  to  update  the  material  stress  state.  The  method  used  here  is  one  developed 
by  Zikry  (1994a)  for  rate-dependent  crystalline  plasticity  formulations.  An  implicit 
finite-element  method  is  used  to  obtain  the  total  deformation  rate  tensor  Dir  To 
overcome  numerical  instabilities  associated  with  stiffness,  a  hybrid  explicit-implicit 
method  is  used  to  obtain  the  plastic  deformation  rate  tensor  Dy.  This  hybrid 
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Parameter 

Value 

°m/°y 

3.0 

Pt/  Ptsat 

0.95 

T/rref 

2.0 

y 

1.0 

a 

2 

b 

1 

for 

0.17 

Table  2.  Properties  of  grains  and  GB  interfacial  regions. 


Grain  bulk 

GB  interfacial  region 

Young’s  modulus  E 

llOGPa 

llOGPa 

Static  yield  stress  <ry 

110  MPa 

330  MPa 

Poisson’s  ratio  v 

0.30 

0.30 

Rate  sensitivity  parameter  m 

0.005 

0.005 

Reference  strain  rate  ]>ref 

0.001s"1 

0.001  s"1 

Critical  strain  rate  ycriticai 

lOV1 

O 

t/3 

Burgers  vector  b 

3.0  x  10"lom 

3.0  x  10“lom 

Reference  stress  interaction 
coefficients  «,■  (;=  1,  12) 

0.50 

0.50 

Initial  immobile  dislocation 
density 

10lom"2 

Varies  as  a  function  of  GB 
orientation,  10lo-1012m"2 

Initial  mobile  dislocation 
density  pf? 

107m"2 

Varies  as  a  function  of  GB 
orientation,  105-107m“2 

numerical  scheme  is  also  used  to  update  the  evolutionary  equations  for  the  mobile 
and  immobile  densities. 


§  7.  Results  and  discussion 

The  multiple-slip  dislocation-density-based  crystal  plasticity  formulation,  the 
specialized  finite-element  algorithm,  the  GB  interfacial  interaction  scheme  and  the 
internal  porosity  formulation  were  applied  to  investigate  the  effects  of  the  interac¬ 
tions  of  multiple-slip  systems  and  mobile  and  immobile  dislocation  densities  with 
GB  regions  to  determine  how  slip  transmission,  impedance,  blockage  and  absorption 
affect  microvoid  nucleation  and  growth  in  a  polycrystalline  aggregate  with  explicit 
void  pairs.  The  material  properties  (table  2)  that  are  used  here  are  representative  of 
polycrystalline  copper  (Zikry  1994b).  A  random-number  generator  developed  by 
Ashmawi  and  Zikry  (2000)  was  used  to  misorient  the  grains  and  GB  regions  by 
the  use  of  the  three  Euler  angles  (p\,d>  and  cp2.  In  this  study,  it  was  assumed  that 
the  GBs  had  random  misorientations  that  did  not  exceed  12°.  GB  regions  were 
assumed  to  have  a  uniform  width  of  0.10  of  the  grain  dimension. 

A  representative  aggregate  size  was  determined  by  modelling  the  response  of 
aggregates  with  different  numbers  of  grains  (20,  30,  50  and  100  grains).  Based  on 
the  convergence  of  the  overall  stress-strain  response,  it  was  determined  that  50 
grains  would  be  sufficiently  representative  of  material  behaviour  for  a  specimen 
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Figure  3.  Polycrystalline  and  random  crystallographic  orientations. 


with  dimensions  of  5  mm  by  10  mm.  This  aggregate  was  subjected  to  an  axial  strain 
rate  of  10_3s_1  by  applying  a  displacement  along  the  [001]  direction  as  shown  in 
figure  3  for  plane  strain  deformations.  Symmetry  boundary  conditions  were  applied 
as  shown  in  figure  3.  Based  on  a  convergence  analysis,  a  minimum  of  1915  four-node 
quadrilateral  elements  were  used  for  the  different  analyses  for  this  study. 

The  initial  mobile  dislocation  density  of  the  grain  bulk  was  chosen  as  107m-2, 
and  the  initial  immobile  dislocation  density  was  chosen  as  10 10  m-2.  In  the  interfacial 
GB  regions,  the  initial  mobile  and  immobile  dislocation  densities  were  varied  as  a 
function  of  the  random  GB  misorientation  (Ashmawi  and  Zikry  2002).  Furthermore, 
it  was  assumed  that  the  initial  GB  static  yield  stress  is  three  times  the  grain  bulk  yield 
stress,  and  that  the  GB  dislocations  are  of  the  same  type  as  the  bulk  lattice.  Using  the 
method  outlined  by  Kameda  and  Zikry  (1996),  the  initial  coefficient  values,  needed 
for  the  evolution  of  the  immobile  and  mobile  densities  given  by  equations  (6)  and  (7) 
were  obtained  as 

ttmmler  —  5.53,  £recov  —  6.67,  ttiminob  —  0.0127,  ttsoiu  —  2.7  X  10  , 

H  ,  (2D 

—  =  3.289  x  103  K. 

K 

Different  arrangements  and  combinations  of  void  pairs  were  used,  such  that  the 
effects  of  void  interactions  and  growth  could  be  realistically  investigated  (figure  4). 
Each  void  was  chosen  as  circular  with  a  radius  of  0.01  mm,  and  each  void  was  placed 
in  the  middle  of  a  host  grain.  Three  void  distributions  were  chosen:  a  horizontal 
distribution  of  two  voids  (case  1),  a  vertical  void  distribution  (case  2)  and  an  inclined 
void  distribution  (case  3).  Each  distribution  had  two  explicit  voids  (figure  4). 
The  voids  were  initially  distributed  0.11mm  apart  for  the  vertical  and  horizontal 
distributions.  For  the  inclined  void  distribution,  the  top  void  was  oriented  at  45° 
from  the  bottom  void,  and  both  voids  were  radially  spaced  apart  by  0.156  mm.  The 
results  of  this  study  will  focus  on  how  the  interrelated  effects  of  saturated  dislocation 
densities  on  different  slip  systems,  dislocation  density  transmission,  impedance  and 
blockage  at  GB  interfaces,  and  void  size  affects  and  controls  failure  initiation  at 
different  physical  scales. 
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Figure  4.  Two-void  arrangements  for:  (a)  case  1,  (b)  case  2  and  (c)  case  3. 


7.1.  Saturated  dislocation  density  effects 

The  saturated  mobile  dislocation  density  for  the  most  active  slip  systems,  as  a 
function  of  the  nominal  strain,  is  shown  in  figure  5.  The  mobile  dislocation  density 
saturation  provides  a  measure  of  how  the  rate  of  inelastic  activity,  such  as  shear 
strain  accumulation  and  porosity,  occurs  on  different  slip  systems.  Mobile  densities 
corresponding  to  the  rotated  void  arrangement  (case  3)  saturated  to  a  maximum 
value  at  an  earlier  strain  than  the  other  two  cases.  Saturation  for  this  case  occurred 
at  a  nominal  strain  of  3.1%  on  the  slip  system  (111 )[10 1],  compared  with  a  nominal 
strain  of  4.5%  for  case  1  (horizontal  voids)  for  the  slip  system  (111 )[1 10],  and  with 
a  nominal  strain  of  4.7%  for  case  2  (vertical  voids)  for  the  slip  system  (11 1  )[1 10]. 

As  can  be  seen  from  the  contours  for  these  mobile  densities  (figure  6),  the 
maximum  mobile  densities  have  accumulated  near  the  peripheries  of  the  voids 
for  case  3.  For  cases  1  and  2,  there  was  no  discernible  mobile  dislocation  density 
accumulation  at  this  nominal  strain.  Furthermore,  as  shown  in  figure  7,  the  accu¬ 
mulated  plastic  shear  was  0.82  at  a  nominal  strain  of  3.1%  for  case  3  compared 
with  0.66  at  a  nominal  strain  of  4.5%  for  case  1  and  0.50  at  a  nominal  strain  of 
4.7%  for  case  2.  These  accumulations  for  case  3  occurred  near  the  void  closest  to  the 
free  surface.  Furthermore,  this  void  has  deformed  into  an  irregular  elliptical  shape. 
In  these  regions  of  accumulated  plastic  slip  and  mobile  dislocation  density,  the 
porosity  has  attained  limiting  values.  As  shown  in  figure  8,  the  porosity  for 
all  three  cases  has  attained  an  approximate  value  of  17%.  Flowever,  the  porosity 
distribution  is  markedly  different  for  cases  1  and  2,  compared  with  case  3.  The 
maximum  porosities  occurred  at  the  top  for  cases  1  and  2.  For  case  3,  the  maximum 
porosity  occurred  in  the  adjacent  regions  near  the  void.  This  indicates  that  microvoid 
nucleation  and  growth  in  terms  of  porosity  /  for  this  case,  can  lead  to  coalescence 
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Figure  5.  Normalized  mobile  dislocation  density  curves  (most  active  slip  systems)  as  a 

function  of  nominal  strain. 

between  the  two  discrete  voids.  Hence,  this  increase  in  porosity  accumulation  is  due 
to  the  void  orientation  along  the  maximum  shear  stress  orientation  of  45°. 
Furthermore,  as  seen  from  these  contours,  for  cases  1  and  2,  the  porosities  have 
accumulated  near  the  free  surface. 

To  ascertain  the  effects  of  the  other  active  slip  systems  on  void  growth, 
coalescence  and  interaction,  the  evolution  of  the  normalized  mobile  density  as  a 
function  of  the  nominal  strain  was  obtained  (figure  9).  As  seen  from  this  figure, 
the  slip  systems  ( 1 1 1 )[  1 10]  and  (11 1)[  101]  saturated  at  much  lower  strains  of  approxi¬ 
mately  3.3%  in  comparison  with  the  slip  systems  (11 1)[101]  and  (11 1)[0 1 1],  which 
saturated  at  a  nominal  strain  of  approximately  7.5%.  This  may  indicate  that  the  slip 
systems  that  saturated  at  lower  strains  may  be  the  operative  slip  systems  that  resulted 
in  initial  pile-ups  at  GB  and  in  the  ligament  region.  These  pile-ups  may  result  in  the 
activation  of  adjacent  slip  systems  that  could  relieve  these  accumulated  densities  by 
transmission  mechanisms.  If  the  pile-ups  are  not  relieved,  then  microvoid  nucleation 
in  these  locations  may  occur.  For  the  slip  systems  (111 )[101]  and  (111 )[01 1]  that 
saturated  at  the  larger  strains,  the  dislocation  densities  were  localized  in  bands, 
in  the  ligament  region,  which  extended  between  the  voids,  where  the  plastic  strains 
have  accumulated  (figures  10  and  11). 

Furthermore,  as  the  global  stress  strain  curve  indicates,  the  stresses  unloaded  at 
approximately  5.1%  (figure  12).  This  is  the  nominal  strain  at  which  the  slip  system 
(11 1)[101]  saturated  (figure  9).  This  further  indicates  that  global  instabilities  are 
directly  linked  to  the  saturation  of  different  slip  systems.  As  shown  in  figure  12, 
maximum  porosities  of  17%  occurred  in  the  regions  of  maximum  mobile  dislocation 
density  near  the  loading  surface  and  adjacent  to  the  point  of  necking. 
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Figure  6.  Normalized  mobile  dislocation  densities  on  the  slip  systems  (a)  (11 1)[1 10]  at  4.5% 
nominal  strain,  ( b )  ( 1 1 1  )[I  TO]  at  4.7%  nominal  strain  and  (c)  (11 1)[10l]  at  3.1% 
nominal  strain. 


Hence,  mobile  dislocation  density  evolution,  saturation  and  interactions  with 
GB  interfaces  play  a  major  role  in  how  dislocation  density  pile-ups  and  transmission 
affect  and  control  void  nucleation,  growth  and  coalescence.  To  emphasize  further 
this  essential  point,  the  transmission,  impedance  and  blockage  of  different  slip 
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Figure  7.  Accumulated  plastic  shear  strains  at  ( a )  4.5%  nominal  strain,  ( b )  4.7%  nominal 

strain  and  (c)  3.1%  nominal  strain. 


systems  at  different  GB  regions  will  be  illustrated  for  the  inclined  void  case  in  the 
next  section. 


7.2.  Grain-boundary  transmission,  impedance  and  blockage 
In  figure  13(a),  a  potential  dislocation  density  system,  corresponding  to  the 
slip  system  (111  )[0 11]  at  5%  nominal  strain  is  shown  for  the  inclined  void  case. 
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Figure  8.  Porosity  distributions  at  (a)  4.5%  nominal  strain,  ( b )  4.7%  nominal  strain  and 

(c)  3.1%  nominal  strain. 


As  seen  from  this  figure,  potential  dislocation  density  activities  were  localized  in  the 
ligament  region,  and  at  the  top  right-hand  side  of  the  aggregate.  These  regions 
are  potential  dislocation  density  regions  that  have  been  identified  by  the  potential 
density  envelope  methodology  outlined  in  §4.  As  the  nominal  strain  increased  to 
8%,  similar  potential  dislocation  density  activity  sites  can  be  activated  (figure  13  (b)). 
Some  of  these  potential  activity  sites  might  be  relieved  by  transmitting  the 
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Nominal  Strain 


Figure  9.  Most  active  saturated  mobile  dislocation  density  systems. 
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Figure  10.  (a)  Normalized  immobile  dislocation  density  for  the  slip  system  (1 1 1)[101]  at 

7.3%  nominal  strain;  ( b )  normalized  mobile  dislocation  density  for  the  slip  system 
(111 )[101]  at  7.3%  nominal  strain. 
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Figure  11.  ( a )  Normalized  immobile  dislocation  density  for  the  slip  system  (1 1 1)[01  l]_at  8% 

nominal  strain;  ( b )  normalized  mobile  dislocation  density  for  the  slip  system  (1 1 1)[01 1] 
at  8%  nominal  strain. 


dislocation  densities  to  the  adjacent  grains,  or  the  dislocation  densities  might  be 
absorbed  by  GBs,  or  the  dislocation  densities  can  accumulate  near  the  GB  regions 
and  evolve  into  pile-ups. 

In  figure  14(a),  the  cross-hatched  area  corresponds  to  a  potential  dislocation 
density  activity  site  on  the  slip  system  (11 1 )[  1 10]  at  8%  nominal  strain.  In  this  case, 
using  the  dislocation-density-GB  interaction  scheme,  96%  of  the  total  dislocation 
density  was  transmitted  through  the  GB,  as  shown  in  figure  14(h),  on  the  slip  system 
(11 1)[  1 10].  From  the  GB,  the  dislocation  density  was  transmitted  by  the  slip  systems 
(111 )[0 1 1],  (111 )[101]  and  (111 )[1 10]  to  the  adjacent  grain  (figures  14(c)  and  (d)). 
Each  of  the  slip  systems  transmitted  one  third  of  the  total  dislocation  density  ema¬ 
nating  from  the  GB.  Again  the  remaining  residual  (4%)  of  the  total  dislocation 
density  in  the  grain  could  evolve  to  a  potential  site  for  dislocation  density  activities. 

These  transmission  and  blockage  interactions  at  the  GBs  are  occurring  during 
the  initiation  and  evolution  of  the  different  deformation  and  failure  modes.  These 
different  interactions  affect  the  evolution  of  porosity  as  the  nominal  strain  changes. 
As  seen  in  figure  15,  the  porosity  distribution  is  affected  by  increases  in  the  nominal 
strain.  At  lower  strains,  the  maximum-porosity  regions  occurred  near  the  void 
peripheries.  However,  as  the  strain  increased  to  8%  nominal,  the  maximum  porosity 
was  concentrated  not  only  in  the  ligament  region  between  the  voids  but  also  in 
the  region  near  the  loading  surface. 


7.3.  Void  size  effects 

To  investigate  the  effects  of  void  size  on  failure  initiation  and  failure,  the  top 
void  size  for  the  inclined  void  case  (case  3)  was  increased  by  100%  to  0.02  mm. 
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Figure  12.  (a)  Global  stress-strain  curve  for  the  inclined  void  case;  ( b )  normalized  mobile 

dislocation  density  for  the  slip  system  (111 )[101]  at  5.1%  nominal  strain;  (c)  porosity 
distribution  at  5.1%  nominal  strain. 


In  figure  16,  the  normalized  mobile  dislocation  density  saturation  curves  are  shown. 
As  the  nominal  strain  evolves  to  8%,  seven  slip  systems  reached  the  saturation  level. 
This  was  in  contrast  with  the  case  with  smaller  voids  (case  3),  where  only  six  systems 
attained  saturation  (figure  9).  These  slip  systems,  for  the  larger-void  case,  saturated 
faster  as  a  function  of  nominal  strain,  in  comparison  with  the  smaller-void  case. 

The  global  stress-strain  response  for  both  this  case  and  case  3  is  shown  in 
figure  17.  The  larger-void-size  case  unloaded  at  a  much  slower  rate  than  case  3. 
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Figure  13.  Potential  dislocation  density  activity  sites  at  (a)  5%  nominal  strain  for  the  slip 
system  (1 1 1)[01 1]  and  ( b )  8%  nominal  strain  for  the  slip  system(l  1 1  )[10 1]. 

This  difference  in  behaviour  is  due  to  the  interaction  of  the  different  dislocation 
density  systems  with  the  GB  interfacial  regions.  In  figure  18,  the  potential  dislocation 
density  activity  sites  corresponding  to  nominal  strain  of  8%  are  shown.  These  figures 
uniformly  indicated  that  the  majority  of  the  transmission  and  pile-up  activities 
occurred  in  the  ligament,  and  near  the  larger  void.  These  interactions  resulted  in  a 
porosity  evolution  pattern  far  different  from  that  for  case  3,  as  shown  in  figure  19. 
For  the  larger-void  case,  at  a  5%  nominal  strain,  these  porosity  distributions  clearly 
indicated  that  void  coalescence  would  occur  along  the  ligament  region  owing  to  void 
growth  and  coalescence.  However,  for  case  3,  the  maximum  porosity  occurred  near 
the  loading  surface  and  at  the  necking  region.  This  difference  in  evolution  is  related 
to  the  presence  of  the  larger  void,  which  resulted  in  different  dislocation-density-GB 
interactions. 

At  a  nominal  strain  of  8%,  although  the  maximum  porosity  for  both  cases 
localized  between  the  voids  (figure  20),  it  was  more  intense  for  the  case  of  larger 
voids.  This  concentration  of  porosity  was  in  both  the  ligament  and  the  necking 
regions.  For  the  smaller  voids  (case  3),  the  porosity  region  near  the  loading  surface 
was  also  localized  in  the  necking  region,  and  it  attained  values  similar  to  those  in  the 
ligament  region. 

In  figure  21,  the  maximum  dislocation  densities  and  the  orientations  of  the 
slip  systems  corresponding  to  specific  slip  systems  were  used  to  delineate  regions 
of  maximum  porosity  at  a  nominal  strain  of  8%.  These  results  emphasized  that  the 
aggregate  with  the  larger  void  will  fail  by  coalescence  in  the  ligament  region  between 
the  discrete  voids.  For  the  aggregate  with  the  smaller  voids,  specimen  necking 
will  most  probably  lead  to  failure.  This  is  consistent  with  the  experimental  observa¬ 
tions  of  Edelson  and  Baldwin  (1962)  and  Bourcier  and  Koss  (1985)  for  void-sheet 
formation. 
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Figure  14.  (a)  Potential  dislocation  density  activity  site  on  the  slip  system  (1 1 1 )[  1 10];  ( b )  96% 

transmission  to  a  GB  on  the  slip  system  (111 )[1 10];  (c)_full  transmission  from  a  GB 
to  a  neighbouring  grain  on  the  slip  systems  (111  )[0 1 1],  (1 1 1)[101],  and  (11 1)[1 10]; 
(d)  distribution  within  the  neighbouring  grain. 


§8.  Summary 

This  paper  provides  detailed  predictive  capabilities  that  have  been  used  to  under¬ 
stand  complex  interrelated  multiscale  physical  mechanisms  that  result  in  void  growth, 
interaction,  coalescence  and  eventual  ductile  failure  in  porous  fee  polycrystalline 
aggregates.  A  multiple-slip  rate-dependent  crystalline  constitutive  formulation 
that  is  coupled  to  the  evolution  of  mobile  and  immobile  dislocation  densities, 
a  new  porosity  formulation  for  microvoid  nucleation  and  growth,  and  specialized 
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Figure  15.  Porosity  distribution  at  (a)  3.1%  nominal  strain,  ( b )  5.1%  nominal  strain  and 

(c)  8.0%  nominal  strain. 


computational  schemes  have  been  developed  to  investigate  the  effects  of  mobile 
and  immobile  dislocation  densities,  dislocation  density  transmission  and  impedance, 
and  GB  misorientation  and  structure  on  void  growth,  interaction  and  coalescence. 
GB  interfacial  kinematic  conditions  have  been  developed  that  account  for  a 
multitude  of  dislocation  density  interactions  with  GBs,  such  as  full  and  partial 
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Figure  16.  Saturation  curves  of  the  most  active  slip  systems  for  the  unequal-sized  void  case. 


Figure  17.  Global  stress-strain  curve. 
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Figure  18.  Potential  dislocation  density  activitysites  at  8%  nominal  strain  for  (a)  (11 1  )[101] 

and  (b)  (lll)[ll0]. 
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Figure  19.  Porosity  distributions  at  5.0%  nominal  strain  for  (a)  the  equal-sized  void  case  and 

(b)  unequal-sized  void  case. 


transmission,  impedance,  blockage  and  absorption.  This  methodology  can  be  used 
to  identify  and  monitor  pile-up  and  transmission  regions  as  deformation  and  failure 
modes  initiate  and  evolve  pertaining  to  the  different  physical  hierarchical  scales 
associated  with  dislocation  density  evolution,  GB  structure  and  orientation,  and 
void  nucleation,  growth  and  interaction. 
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(a)  (b) 

Figure  20.  Porosity  distributions  at  8.0%  nominal  strain  for  (a)  the  equal-sized  void  case 

and  ( b )  unequal-sized  void  case. 


Figure  21.  Failure  traces  at  8.0%  nominal  strain  for  (a)  the  equal-sized  void  case  and 

( b )  unequal-sized  void  case. 


For  aggregates  with  different  void  pair  arrangements,  it  was  shown  that  mobile 
dislocation  density  evolution  and  its  saturation  on  different  slip  systems  and  its 
interaction  with  GB  interfaces  play  a  major  role  in  how  dislocation  density  pile- 
ups  and  transmission  affect  and  control  void  nucleation,  growth  and  coalescence. 
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When  the  mobile  dislocation  density  saturates  on  a  specific  slip  system,  this  provides 
sites  for  localized  interactions  between  explicit  voids  and  substantial  increases 
and  localization  of  porosity  and  accumulated  plastic  strain.  However,  the  mobile 
dislocation  density  saturation  may  also  lead  to  GB  blockage  or  absorption  within 
the  GB  region.  This  can  occur  because  there  may  not  be  preferentially  oriented 
slip  systems  in  adjoining  grains  to  relieve  a  build-up  of  the  accumulated  dislocation 
density  activities  in  the  GB  regions. 

Analyses  were  also  conducted  to  determine  how  changes  in  void  pair  sizes  and 
spacing  affected  void-to-void  interactions,  and  how  this  was  related  to  dislocation 
density  interactions  with  GB  interfaces.  It  was  shown  that,  as  the  void  pair  sizes 
were  increased,  the  porosity,  the  stresses  and  the  accumulated  plastic  strains  were 
localized  in  the  ligament  region  between  the  voids.  This  occurred  because  there 
was  a  significant  number  of  regions  for  potential  dislocation  density  activities 
as  indicated  by  the  saturation  of  the  mobile  dislocation  densities  of  different  slip 
systems  in  the  ligament  region.  As  the  void  size  was  decreased,  there  were  fewer 
saturated  slip  systems  and  activities  in  the  ligament  region.  There  were  also  higher 
local  stresses  associated  with  smaller  voids.  The  stresses  are  higher  for  the  smaller 
void  sizes,  because  of  the  build-up  of  dislocation  densities  at  the  GBs  and  around 
the  void  peripheries,  which  eventually  leads  to  the  global  unloading  of  the  stress. 
Therefore,  there  are  two  different  types  of  failure  associated  with  void  pair  size. 
As  the  void  size  is  increased  for  different  void  pairs,  porosity  and  interaction  are 
localized  in  the  ligament  region.  As  void  size  is  decreased,  failure  is  due  to  the 
global  unloading  of  the  stress  as  a  result  of  build-up  of  local  stresses  adjacent  to 
the  void  peripheries. 
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